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>^ Abstract 

The calculation of Franck- Condon factors between different one-dimensional Morse 
potential eigenstates using a formula derived from the Wigner function is dis- 
cussed. Our numerical calculations using a very simple program written in Math- 
ematica is compared with other calculations. We show that our results have a 
similar accuracy as the calculations performed with more sophisticated methods. 
We discuss the extension of our method to include non-Condon effects in the 
calculation. 
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1 Introduction 



The study of complex features of molecular vibrational spectra at high overtones and combi- 
nations has become possible by recent advances in optical laser spectroscopy Rovibrational 
level structure of small polyatomic molecules has been elucidated in a formerly inaccessible 
energy range. Transition intensities between these levels can reveal fine aspects of the struc- 
ture and dynamical properties of these systems. Theoretically, the evaluation of molecular 
vibrational wave functions has relied upon a number of assumptions and approximations, 
of which the use of a harmonic normal-mode basis is the best known procedure. At higher 
energies, however, normal-mode vibrational assignments of individual states becomes more 
difficult due to anharmonic effects. In this context, Morse potentials |]J have been proposed 
to more accurately model molecular systems, as they often constitute good approximations 
to diatomic molecular potentials derived ab initio M. Even if the general solutions to the 
corresponding Schrodinger equation cannot be written in analytical form, the case of zero 
angular momentum can be closely approximated by the one- dimensional analytical solutions 
@. The analytic form of these wave functions has been used not only for the simple diatomic 
case, but for the analysis of vibrational structure in polyatomic molecules, within a local- 
mode framework (see eg. @, 0, ^, 0)- In addition, the connection of the one- dimensional 
Morse, Poeschl- Teller and other exactly solvable potentials with the SU(2) group has been 
exploited in the formulation of algebraic models which incorporate the anharmonic character 
of the local vibrations from the outset (see || and references therein). 

In recent studies, the Lie-algebraic approach has been also applied to the evaluation of 
Franck-Condon intensities ||, [L0| (FC) in small polyatomic molecules 0, |^, 0. These studies 
accurately reproduce the emission spectra of the S2O C X A — X X A! (vr* — 7r) experimental in- 
tensities ||, requiring for their calculation a combination of two procedures: a) an accurate fit 
of the vibrational energy levels in both electronic states using a local-mode algebraic Hamil- 
tonian, and b) an efficient and accurate methodology for the evaluation of one-dimensional 
Morse wave function overlaps. One possibility to calculate overlaps of Morse wave functions 
is to use an approximation based on an energy-dependent harmonic approximation where the 
widths and displacement of the harmonic potentials are adjusted in order to mimic the be- 
havior of the Morse potentials at the corresponding energy levels [|J . Although this method 
gives explicit formulas for the overlap integrals, and consequently a good accuracy for the 
low lying vibrational states, it departs from the exact results at higher energies. Other (nu- 
merical) methods, like the use of the Configuration Localized States (CLS) formalism and 
the use of Gaussian quadratures [[□], [TJ|] have been proposed as alternative ways to evaluate 
the overlaps of two Morse functions. 

In view of future applications to the evaluation of multidimensional Franck-Condon factors, 
a simple and reliable method of calculation of one-dimensional Morse functions overlaps has 
to be devised. In this paper we propose a new method of evaluation based on our recent work 



on the Wigner function of the Morse potential eigenstates [|13|]. Our main result contains 
analytical sums of simple one-dimensional integrals which are not tabulated and which must 
be evaluated with standard numerical techniques. We make no approximations but must 
rely on numerical evaluations for part of our calculations. Fortunately, present computa- 
tional resources (like Mathematica) can handle these calculations with very high accuracy in 
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a simple and efficient way. We show that a short program written in Mathematica (included 
in an appendix) gives accurate results for several molecular examples. To demonstrate the 
reliability of the method we compare our results with other calculations and provide esti- 
mates for its range of applicability. In the last section we show that the procedure can be 
extended, with no loss of accuracy, in order to include non-Condon effects. 



2 The Morse Potential Eigenstates 

In what follows we briefly review the properties of the Morse potential eigenstates and the 
use of the corresponding Wigner function to evaluate their overlap. We then express our 
result in terms of the integral J a (eq. (|22|)), for which we find a symmetry property and 
a recurrence relation, which can be used to verify the accuracy of its numerical evaluation. 
After discussing some particular analytic solutions, we apply our method to several examples 
and compare our results with more sofisticated techniques [|[ |ll|, [Hfl . 

The Morse potential is characterized by two parameters: its depth T> (bonding energy), and 
the inverse of its range (3 > (restitution constant in the harmonic approximation). In the 
one-dimensional version it is given by 

V{x) = d\(1 - e- 13 *) 2 - l\ , (-oo<a;<oo), (1) 
where x is the position with respect to the equilibrium point. 

It is well known that the discrete spectrum of the Morse potential is described by the formula 
where 

. V^D 1 

while v is the number of anharmonic phonons u = 0,1, . . . \J\, with being the largest 
integer not exceeding j, and \i the reduced mass of the molecule. 

If j takes an integer or half integer value, it determines the irreducible representation Qi of 
the SU(2) algebra which is the dynamical algebra of the Morse potential. In general there 
are \J\ + 1 bound states, except when j is an integer, in which case one state has zero energy. 
The eigenfunctions corresponding to the eigenvalue (0) are given by 

C(0=^e-«/ 2 e^ '^(0, (4) 
where a new variable (the Morse coordinate) 

£ = (2j + l)e-^, (5) 
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has been introduced. The wave function (f|) includes the associated Laguerre polynomials 
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normalized by the condition 

poo 

d*|<(*)| 2 =l, (7) 



-oo 

when the normalization factor is taken in the form 



{Kf _ /M2(j-„) 



r(2j - v + 1) 

The eigenvalues (|]) can be written as 

E u = —D + hw e (v + 1/2) - hw eX {v + 1/2) 2 , (9) 

where 



/2D 

hw e = f3hJ— , (10) 

h 2 [3 2 

k, e x = — . (11) 

These relations relate the parameters of the Morse potential with the standard spectroscopic 



parameters ui e and uo e \ ||16||. In particular we have 



2 J + 1 = ^ = 4^ , (12) 



and 



> \\^[LUJ e x ■ (13) 



3 Franck-Condon factors for the Morse Potential 



Molecular Franck-Condon factors describe, in first approximation, the dependence on the 
vibrational wave functions of transition intensities between vibrational states in different 
electronic configurations. Assuming that the electronic wave functions are independent of 
the vibrational states (Condon approximation [|Kj]), the Franck-Condon factors are given by 
the square of the overlap integral between initial and final wave functions in the transition 

*i,2 = |/i, 2 | 2 (14) 
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where 



fi,2 = (^(x - Ri)\il%„(x - R 2 )) , (15) 

and R\ and R2 are the equilibrium distances of the nuclei in the initial and final electronic 
states, respectively. Relative transition intensities, in the Condon approximation, are simply 
related to the Franck-Condon factors flip) by || 



h,2 oc z^ 4 |/i,2 1 2 , 

where the frequency dependent factor v A (y is essentially the difference between vibrational 
energies of the initial and final states in the transition), is introduced for calibration with 
experimental data. 

In a higher approximation (Non-Condon), the electric dipole operator for transitions is ex- 
pressed as a series expansion on the interatomic coordinates |T7|] or, equivalently, as a series 
expansion in one of the Morse coordinates £,\(x) or ^(x) Q. Our method can be adapted to 
include such effects in a straightforward manner (see discussion below). 
A simple way to calculate the FC factors arises from the Wigner function of two Morse wave 



functions, which we define in the form JT3 



1 r°° 1 1 

W K\^^U^p) = —J drifts {x - -r) e-^^Jx + -r) . (16) 

Making the change of variable z = \r and taking into account that 



3l"l , 

where ipj^i ( x ) * s the wave fuction corresponding to the mirror image of the Morse potential 
([!]) whose normalization is clearly also given by (|§|). 

We then find the expression 

W{^l%^\x,p)=~ I dz^ Ui \x-z)e- 2 ^ h ^l V2 (z + x). (18) 

If we evaluate the integral at (x = |(i?2 — R 2 )^P = 0) we get precisely the overlap integral 
between and Vv^' and thus 

/i,2 = -nhW(i;-%,^l U2 \x = I(i2a - R 2 ),p = 0) . (19) 

In Appendix A we calculate the Wigner function corresponding to two different Morse po- 
tential states (H), leading to the following expression for the Franck-Condon factor: 



5 These coordinates are the natural coordinates appearing in the eigenfunctions of the Morse potential 
and are defined as £i(x) = (2ji + 1) exp (—fliXi) and £2(2;) — (2j2 + 1) exp (—P2X2) (see Eq. @) 
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Ji ~ \l T (2j 1 -u 1 + l)T(2j 2 -u 2 + l) yi V2 



X 

^1/2=0 



x J a (ji - vx +£ u j 2 - u 2 +4,2/1,2/2) , (20) 



where 



and 



yi = (2j 1 + l)e-^ R2 - R iV 2 , 

y 2 = (2j 2 + l)e +l3 ^- R ^ 2 , (21) 

f°° ,,71 --(Vl^ + V2« a ) 

Uh,k 2 ,y u y 2 )= duu k ^- l e 2 KUl \ (22) 

Jo 

with a = 1^, fei = ji - ivi + 4 and fc 2 = h — v i + h- 

Overlap integrals of the type ( p^) were previously encountered in ref.(18|], where it is stated 
that they "cannot be evaluated exactly in general (for a 7^ 1)" and "it is necessary to resort 
to a valid approximation" . In addition, it was stated that the alternating sums involved in 
( P0|) may give unstable results as a consequence of the numerical error accumulation when 
the sums are performed over very large numbers, a situation that can arise for realistic 
calculations involving large values of j\ or j 2 . 

As mentioned in the introduction, several approaches have been proposed to find FC fac- 
tors (|20|) in analytic form, including the use of different versions of the harmonic oscillator 
approximation of the Morse wave functions j|, |^, [19], ^0], the expansion of the integral 



(22) into series containing polygamma functions the use of a Laguerre quadrature [11 



and others. We now show that the integral fl2"2]), although not tabulated, is an analytical 



well-behaved function which satisfies relations which can be used for its reliable evaluation 
using a personal computer. In particular we have found that the program Mathematica can 
be applied to calculate those integrals with high accuracy in an efficient way. 



We start our study by stating the following properties of the integral J a : 
1. Symmetry property: 

%{ki, h, Vi, 2/2) = - 3i/a{k 2 , h, 2/2, 2/i) , (23) 
a 

which follows from fl22| ) and the change of variables u' = u a . 

As an application of (|23|) we can check the accuracy of the integration appearing in 
the transition A x E+(u = 0) - X l Y>+(v = 0) of the 7 Li 2 molecule (Table f§. In 
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this case the difference of both calculations 3 a — - 3\i a is —1.71551 x 10~ 67 where 
3 a pa - Ji/ a pa 1.14752 . . . x 10~ 53 . This shows the typical precision in the computations 
of integrals (0). In order to get such an accuracy we have used a scaling method for 
the integral. This procedure is discussed below. 

2. Recurrence relation: 

3 a (h,k 2 ,y 1 ,y 2 ) = —[2(^i + ak 2 - 1)3 a (h - l,k 2 ,y%,y 2 ) 

-ay 2 3 a (k 1 - l,k 2 + l,y!,y 2 ) + 5 aA ] , (24) 

where a = k\ + ak 2 . This property follows by carrying out the integration by parts. 
Note that a = 1 can only occur for the case a = 1 (f3\ = /3 2 ). 



Regarding the application of the three-term recurrence relation fl24|) , we observe that it 
can be used to find one such integral if the other two included in this relation are known. 
It is therefore possible to find a full set of integrals (|22|) with all possible values of ki = 
1, . . . Ni, k 2 = 1, . . . N 2 if the integrals 3 a (k 1 — 1 . . . N x + 1, k 2 — 1, y 1: y 2 ) and 3 a (ki = 1, k 2 = 



1 . . . N 2 + 1, yi, y 2 ) are known for given a, y\, and y 2 . The relation (|24D can be used in two 
ways: to check the accuracy of the numerical integrations and to calculate the values of the 
other integrals. 



As example we use the recurrence relation (|24|) to check the accuracy of the integrations 
involved in the transition A 1 S+(z/ = 6) — X l Y^(v = 6) of 7 Li 2 molecule (see Table ^). 
Let us consider the integral 3 a (Eq. |22|) with (£i = 0,£ 2 = 1) and compute it by numerical 
integration and with the recurrence relation (p4|). For this we need to compute integrals 
with {t\ = 0,£ 2 = 0) and {t\ = 0,£ 2 = 1). The absolute difference of both calculations is 
3.17601 x 10~ 67 , being 3 a ~ 2.07527. . . x 10~ 53 , in agreement with the estimated precision 
obtained with the symmetry relation. 

Relations (^) and (|24]) led to a simple method to verify the accuracy of the numerical 
integration. We are currently studying other properties of the function (p2[), which is a gen- 
eralization of the T function, to which it reduces in the case a = 1. 

In realistic cases, quite often the overlaps must be calculated for states with rather large val- 
ues of the parameters j% or j 2 (or equivalently for large bonding energies, see Eq. (^). In this 
cases the numerical integration of (^) requires a special treatment in order to achieve high 
accuracy. Although the integrand is a well behaved function of x having a single maximum 
and going to zero at x — > and x — > oo, the resulting integral, being a very small number, 
could be difficult to evaluate due to the limited precision of the numerical calculations. One 
way to deal with such complications is to perform a change of variables in order to scale the 
integrand in such a way that the integrand at the maximum acquires a unit value at x — 1. 
In all cases studied, we have found that the best procedure to evaluate (|22|) is by carrying out 
a scaled integration. Our tests showed that the relative differences of the integrals calculated 
using formula ( |2"2"| ) and the symmetry relation (^) are very small (< 10 -13 ) when the scaled 
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integrations are used in both cases. 



As our first example we have calculated Franck-Condon factors for the S — S mode in the 
5*20 molecule (Table |I]). We found a maximal relative difference ~ 3.14 x 10~ 14 between the 
scaled integrals ( p2]) and the scaled integration using the symmetry relation (|23|). Despite 
the fact that the numerical integration can be calculated with high accuracy, it is not always 
enough to compensate for the loss of accuracy due to the alternating sums involved in (ED) . 
This becomes more apparent in the cases in which we have a significant number of terms 
in the summations. We can see this effect in the results of Table |T] (e.g. for v% = 7 and 
z/g = 6). The use of logarithms to evaluate the products (and divisions) of the very large 
numbers appearing in (|20|) is usually a good strategy to avoid some problems. However we 
have also calculated the FC factors using this strategy and found no significant improvement 
in the results. 

In Table [L] we compare our results with calculations using the method of quadratures pro- 



posed in [O, 12 1, and with the algebraic formulation using modified harmonic functions of 



ref. H. Our results are in excellent agreement with the results of ref. |TT[ except for the 



corner in the table close to = 7 and Vq = 6 where the calculations exhibit some deviations 
due to the above mentioned loss of accuracy in the sums. 

The approach using modified harmonic functions is in this case expected to be a good ap- 
proximation (at least for the overlaps between the lowest-lying states) since we have large 
values of the parameters ji,j 2 (large depths Di,^) • However, there are significant devia- 
tions with respect to our results even for this extreme case. 



For the sake of completeness we now discuss some exactly solvable examples of the FC for- 
mula @g): 

1. Pi = Pi (a = 1). In this case we have fTT|j: 



J =i(ji + t\,h -"2 + 4, 2/i, 2/2) 



duexp 



-^(2/1 + 2/2)^ 



u 



a-l 



T(a) 



2/i + 2/2 



with a = ki + k 2 = ji — v x + ii + j 2 — v 2 + £2, and thus 

2, 



Jl,2 



VlW{jl ~ ^l)(j2 - V 2 ) 



r(2 jl - v x + i)r(2j 2 -i/2 + i 



x E 

^1,^2=0 



2ji - v\ \ ( 2j 2 - u 2 



xyi £l V2 h 



l x \t 2 \ V vx-l\ ) \ U2-I2 



T(a) . 



2/i +2/2 



(25) 



(26) 
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If we define: 



c 



2j 2 



I 



-P(R 2 -Ri) 



m 
yi 



(27) 



2ji + 1 

we obtain the same expression as Eq. (4.5) in Ref. The results for the Franck- Condon 
factors obtained trough the numerical integration of ( [22]) and with the analytic formula ( [ZED 
are of comparable accuracy. This gives another indication that the main source for the loss 
of accuracy is due to the alternating sums appearing in (pOp, (|26|) . To check the accuracy of 
the numerical integration (|22"D, we have calculated its relative difference with the analytical 
result (|^) for a broad range of values (e.g. for 2/1 = 2/2 = 100 and a = k± + k 2 — 1 . . . 200 
) and found a maximal relative difference of A(J a ) < 5 x 10~ 13 . However we should note 
that the analytical result for the integration ( ^5[) can also lead to uncertainties when it is 
numerically evaluated, specially for large values of a and/or 2/1,2/2- 



2. p 2 — 2/3i, (a = 2). In this case we have 
Ja=2(ji -vi + £1,32 -v 2 +h,yi, 2/2) 



du exp 



a 



2/2 



2 r, 



a) exp 



162/2 



D_ a being the cylindrical functions with a — j\ — V\ 



Ayiu + y 2 u 1 



Vi 



u 



0-1 



(28) 



2^ 

'1 — 2j2 — 2z/2 — 2^2, and thus 



ra=2 
J 1.2 



) , 2z/ 1 !z/ 2 !(j 1 -z/ 1 )(j 2 -^) ^ j2 _, 2 

r(2ji-^ + i)r(2j 2 -i^ + i) w y2 



€lr«2=0 



?i!£ 2 ! V ^1 



^2 - « 2 



a 



2//V 2 



(2/2) 2 n 



a J exp 



162/ 5 



2/1 



2^ 



(29) 



3. P2 — 20i, a — \ 
In this case we have 



J a =i (Ji - ^1 + ^1, J2 - ^2 + £2, 2/i, 2/2) 



exp 



-^(yiu + y 2 Vu) 



u 



a-1 



= 2J2O2 - ^2 +4,jl - ^1 + ^1,2/1,2/2) 

where we have used the symmetry property (p3|) , and therefore 



(30) 



/I. 



2z/i!z/ 2 !(ji - ^i)(j 2 - ^2) 
T(2j 1 - v x + l)r(2j 2 - ^2 + 1 



2v ^w D Am • (3i) 



with a = 2ji - 2ui + 2l x + j 2 -u 2 + i 2 . 

These particular cases can also be used to check the accuracy of the general result (5(1 122|) . 



4 Examples 

Continuing with our comparison with other approaches, we have considered (i) the case of 
two identical (but displaced) Morse potentials and (ii) the case of two different Morse po- 
tentials. We have taken these examples from ref. [O]. Our results are identical to those 



obtained in |TT] where an analytical approximation for the Morse functions overlaps based 
on the use an integration by quadratures was used Q (see Tables (fj) and (|3|)). Perhaps this 
is not surprising since the corresponding values of ji,j 2 are small. The accuracy of the nu- 
merical calculations is high given that the effects of the numerical instabilities coming from 



the alternating sums are relatively small. The relevant point, as observed in is that the 
deviations in the results obtained from harmonic functions relative to those using the Morse 
functions are larger when the depths D (and correspondingly the j values) are small. 

As a next set of examples we have considered the vibrational transitions studied by Ley-Koo 
and collaborators Hl4| . In ref. [J4] Franck-Condon factors of diatomic molecules were calcu- 



lated by means of 3-dimensional Morse functions via the method of confinement in a box. 
The solutions of the radial Schrodinger equation were found by diagonalization of the Hamil- 
tonian in a large basis (N ~ 600) of free-particle eigenfunctions with appropriate boundary 



conditions. We have borrowed some tables from ref. |L4j to facilitate the comparisons. Our 
results for the Franck-Condon factors are in excellent agreement with those of Ley-Koo and 
collaborators (see Tables f|, [5], |, 0). This shows on the one hand that in spite of the fact 
that we use 1-dimensional Morse functions, the contributions coming from the unphysical 
region (— oo < x < 0) are entirely negligible, and on the other, that our simple procedure 
achieves results with a similar accuracy as those arising from the numerically intensive and 
sophisticated methods. 

The specific transitions which are included for comparisons are: 
1. A 1 ^ - X X S+ in 7 Li 2 molecule (Table §. 

6 In Ref. [JO] it was suggested to calculate the Franck-Condon matrix elements using the Laguerre 
quadrature in order to avoid the calculation of the alternative sum in (p0|). However this approach may 
present a similar disadvantage since the calculation involves a sum over all positive and negative nodes of a 
polynomial. 
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2. B X U+ - X X S+ in 7 Li 2 molecule (Table §. 

3. A 2 U 2 - X 2 S+ in CN molecule (Table |) . 

4. B 3 U g - A 3 £+ in N 2 molecule (Table 0). 

Tables |], |5|, ||, [7] contain results of calculations arising from analyses of experimental data 



22, 23, M 



], simple harmonic oscillator approximations |2j 



imations 20 



|, anharmonic approx- 
29], graphical integration pp |, asymptotic expansions |^TJ and by the method 

We emphasize that we are concerned with 



of confinement of Morse potentials in a box |14 



an accurate evaluation of the overlap of two Morse functions and thus comparisons have to 
be done with similar calculations, in particular with the results obtained by Ley-Koo and 
coauthors f|T3 . 



5 Non-Condon Effects 

Our method can be generalized to include non-Condon effects by taking into account more 
terms in the expansion of the electric dipole operator. Assuming a power expansion in terms 
of one of the Morse coordinates, say £i [|3~2 



r?=0 



(32) 



and following a similar calculation to that leading to (|20|), we obtain for the dipole transition 



(j2,^2|T|ji,z/i) = / M (x)7(£x(x))il> hn (x)dx (33) 



]2 V 2 

— OO 



2ji - v\ \( 2j 2 - V2 



Vi 1 V2 i2 x ^2a v - h a (k 1 +r],k 2 ,y 1 ,y 2 )> 

?7=0 ^ ' 



with a = k\ = ji — v\ + i\ and k 2 = j 2 — v 2 + £ 2 . The term r] = corresponds to the 
FC part. The only difference with ( f2~0D is that instead of a single integral J a (&i, k 2 , y\, y 2 ) we 



have a sum of integrals 3 a (ki + rj, k 2l y\ 1 y 2 \ each of them multiplied by the corresponding 
coefficient in the expansion ([32]). That means that we only have to calculate some more 
integrals and multiply them by some coefficients which should be chosen in a suitable way. 

In Figure 1 we have plotted the calculated emission intensities for the transition C l A' — X l A' 
of the S 2 molecule in the SS mode using (a) a zero order approximation (FC), and (b) 
including the first correction in (|33|) for some values of the coefficient a\. The emission 
intensities were calculated for the transitions (0 — > v) and (1 — > v). We can see that the 
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effect of the first non-Condon term reduces, as expected, the transition intensities to the 
higher levels (large v) being more important for the transitions (1 — > v). Although we 
have not attempted to perform a detailed analysis of this transitions (since in any case a 
full polyatomic treatment is necessary), it is clear that the present approach can be easily 
adapted to include non-Condon effects without major changes in the method. 



6 Summary and Conclusions 

In this paper we have calculated Franck-Condon factors by means of a simple formula ( ^0[) 
(derived from the Wigner function) which can be adapted in a straightforward manner to 
include non-Condon effects when the transition operator is written as a series expansion in 
powers of one of the Morse coordinates (|33|). The formula obtained involves well behaved 



one-dimensional integrations ( |22|) which are not tabulated, but for which a symmetry prop- 
erty (^3|) and a three term recurrence relation (|24]) have been derived to check the accuracy 
of their numerical evaluation. We found that the best strategy to treat the numerical inte- 
grations is by means of a scaling method in which the integral is transformed in such a way 
that the integrand acquires a unit value at the maximum shifted at x — 1. This approach 
is particularly useful when the transitions involve large bonding energies (or large values of 
the parameters ji,]?)- 



Since the present approach is oriented to possible applications to the calculation of multi- 
dimensional Franck-Condon factors, we have confined our study to the calculation of one- 
dimensional Morse function overlaps by using a simple program written in Mathematica (see 
appendix B for a listing of the Mathematica program used to calculate the FC factors). 



As far as we only require transitions between the lowest-lying states in each Morse potential 
(up to V\ « u 2 ~ 10) our results have a similar accuracy as those obtained through more 
sofisticated methods. Although in this work we have included only a few examples, we were 
able to reproduce all results in [14| with high precision. Only for cases involving long summa- 
tions we observe some isolated small deviations. We emphasize that in the above transitions 
we are dealing with rather large values of the bonding energies in the Morse potentials (i.e. 
large summations in eq. (|20"1) ), which suggests that our procedure can be reliably applied to a 
wide variety of physical examples. On the basis of these results, it is our impression that the 
inaccuracies arising from the alternating sums in (|20|) have been overestimated. Such sums 
can be calculated by Mathematica with no significant loss of precision for a wide range of 
problems of physical interest. However, for transitions involving higher-states a more elab- 
orate programming of our formula is required in order to avoid the loss of accuracy. This is 
beyond the objective of the present work and will be studied elsewhere. 



In conclusion, Franck-Condon factors for the Morse potential and corrections including non- 
Condon effects can be calculated with high accuracy by means of a very simple program 
suggesting that this approach is very convenient, since it can be implemented in any per- 
sonal computer. 



11 



A Derivation of Franck- Condon overlaps for Morse Po- 
tential eigenfunctions 



Let's consider the Wigner function of two different Morse wavefunctions if)j V (x) and ip 



p2 



2nh 



■■'-■■"^ / ,/rox,)«| -I(2j 1 + l)e PH * 2 n } 



(2jt + l)e 



ji-i'i 



(2j! + l) e - A ( x -3 r )' 



-ipr/h 



exp --(2j 2 + l)e 



-/3 2 (x + -r) 



(2j 2 + l)e 



-/5 2 (x + -r) 



.72 -^2 



'VI 



(2j 2 + l) e -A(*+|r) 



defining 



2/! = (2j! + l)e- fta \ y 2 = (2j 2 + l)e 



■fax 



and making the change of variables 



u = e 



e P 2 r/2 = uWl du= f^ udr 



we have 
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_ju^_h^2 (ii - y\ ) (j2 - ^2 ) 
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00 11 2 ' 

-oo 2 2 

"2. 



L^-^( yi n)L^-^(y 2 n-ft 

then, expanding the associated Laguerre polynomials in their finite series (|6]), we arrive 
finally to 



g (_1)^2 



^1^2=0 



Wit 



V-Z ) ( 2 £-Z i 



1 1 « /« C/i — ^i — 75-U2 - ^2) - + 1\ - -rh - 1, 
duexp{ - -y x u - -y 2 «~ ft/ft }^ ^ ^ 



from which we obtain the Franck-Condon factors 
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B The Mathematica program for the calculation of the 
Franck- Condon factors with wave functions of the 
Morse potential. 

FCFact [ j 1_ , nul_ , bet al_ , Rl_ , j 2_ , nu2_ , beta2_ , R2_ , f lag_] : = 
Module [ 

{yl,y2,a}, 

(* *) 

(* This program computes the Franck-Condon factor of two Morse functions *) 
(* Usage: *) 
(* FCFact [jl,nul,betal,Rl, j 2, nu2,beta2,R2, flag] *) 
(* flag=0 computes the FC factor *) 
(* flag=l computes the square of the FC factor *) 

(* *) 

(* Definition of Local Variables *) 

Acc=12; 

AccN=100; 

a= SetAccuracy [beta2/betal , AccN] ; 

yl=SetAccuracy[(2*jl+l)*Exp[-betal*(R2-Rl)/2] ,AccN] ; 
y2=SetAccuracy[(2*j2+l)*Exp[-beta2*(Rl-R2)/2] ,AccN] ; 
Fcsum=0 ; 

f acOa=SetAccuracy [2*Sqrt [a*FullSimplif y [nul ! * ( j 1-nul) *nu2 ! * ( j2-nu2) / 

(Gamma[2*jl-nul+l]*Gamma[2*j2-nu2+l])]] ,AccN] ; 
facOb=SetAccuracy[yl~(j 1-nul) *y2~(j2-nu2) ,AccN] ; 
f acO=f acOa*f acOb; 
Do [ (* Do loops for summations *) 
f acla=Set Accuracy [(2* j 1+1) " (11) * (2*j2+l) ~ (12) , AccN] ; 

f aclb=SetAccuracy [Exp [-betal* (R2-R1) * (11) /2] *Exp [-beta2* (R1-R2) * (12) /2] 
,AccN] ; 

f acl=SetAccuracy [f acla*f aclb , AccN] ; 

fac2=SetAccuracy[(-l)-(ll+12)/(ll!*12!) ,AccN] ; 

f ac3=SetAccuracy [FullSimplif y [Binomial [2* j 1-nul ,nul-ll] * 

Binomial [2* j2-nu2 ,nu2-12] ] ,AccN] ; 
Fcaux=SetAccuracy [f acO *f acl*f ac2*f ac3,AccN] ; 

(* =========================================================== *) 

(* Integration Section. Integral la. Calculation with scaling *) 
( * =========================================================== * ) 

Lambda=SetAccuracy [jl-nul+11 + a*(j2-nu2+12 ),AccN]; 
(* Integrand *) 

integr=x~ (Lambda- 1 ) *Exp [- (yl*x+y2*x~a) /2] ; 
finteg[xx_] :=xx"(Lambda-l)*Exp[-(yl*xx+y2*xx~a)/2] ; 
(* Looking for position of maximum *) 

aux=FindMinimum [-integr , {x, 1,0, 1000} ,PrecisionGoal->10] ; 
peak=aux [ [2] ] [[1]] [[2]] ; 

xatmax =Set Accuracy [peak, AccN] ; (* Position of maximum *) 

(* Integrand value at maximum *) 
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f atmax=SetAccuracy [f integ [peak] ,AccN] ; 
finteg2=Set Accuracy [(1/fatmax) *(X*xatmax) " (Lambda- 1) * 
Exp[-(yl*(X*xatmax)+y2*(X*xatmax)~a)/2] ,AccN] ; 
res=Set Accuracy [f atmax*xatmax*NIntegrate [f integ2 , {X, , Infinity} , 
PrecisionGoal->Acc , AccuracyGoal->Acc , 
Method->GaussKronrod] ,AccN] ; 
(* Adding terms *) 

Fcsum=Fcsum+SetAccuracy [Fcaux*res , AccN] 

,{ll,0,nul},{12,0,nu2}] ; 
(* Choosing Output *) 
lf[flag==0, FCf actor=Fcsum] ; 
If[flag==l, FCfactor=FcsunT2] ; 
N[FCfactor,15] 
] 
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Figure 1: Emission intensities in the S2O C X A — X X A transitions (a) and (b) I\-* v . 
The calculations in the Condon approximation using Eq.(p0"D (FC) as well as the calculations 
with inclusion of the first non-Condon correction as given by Eq. (|33|) (NC) for a\ = 0.01 
and a\ = 1 are displayed. Data from experimental analysis indicated by crosses is from 
ref.||. The parameters for this calculations are ji = 128, /5i = 1.6, j'2 = 79, $2 = 1.8 and 
A = i? 2 - -Ri = 0.26 
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Table 1: Franck-Condon factors (^) for v%_ = ... 7 and vq = . . . 6 for the S — S stretching 
degree of freedom in the S2O molecule. Parameters for overlaps are ji = 128, /3i = 1.6, 72 = 
79, fa = 1-8 and A = R 2 — R\ = 0.283. Successive entries correspond to the calculations 
using formula (|2~0|) with a scaled integation (first line), using the Configuration Localized 
States method of ref.jll, 12 and using the algebraic formula of ref.[Q. 



fl\^2 





1 


2 


3 


4 





0.42687 


-0.34228 


0.27839 


-0.21679 


0.14889 




0.38163 


-0.52971 


0.51990 


-0.41663 


0.28915 




0.38163 


-0.45147 


0.49396 


-0.52002 


0.50695 


1 


0.79048 


-0.12775 


0.04693 


-0.00110 


- 0.01166 




0.52971 


-0.35362 


-0.02750 


0.32220 


-0.43193 




0.62781 


-0.28358 


0.05443 


0.16330 


- 0.38378 


2 


0.41678 


0.74799 


-0.13014 


0.13408 


-0.08036 




0.51990 


0.02750 


-0.38061 


0.28255 


0.03174 




0.23717 


0.51298 


-0.40591 


0.24204 


0.01297 


3 


-0.12667 


0.44613 


0.77692 


0.04596 


0.09256 




0.41663 


0.32220 


-0.8255 


-0.15417 


0.35170 




-0.09261 


0.32338 


0.41678 


-0.40496 


0.29047 


4 


0.02807 


-0.24277 


0.21095 


0.81696 


0.36737 




0.28915 


0.43193 


0.03174 


-0.35170 


0.08991 




0.05862 


-0.19434 


0.36082 


0.38247 


-0.39483 



Table 2: Franck-Condon overlap integrals fi.2 between states in identical Morse potentials 
displaced one with respect to the other. The Morse parameters are j\ = j'2 = 5, /3i = fii — 
0.90 A , Ri = 2.67A and R^ = 3.60A. The first line correspond to the results of formula 
(Iffi) using numerical integrations fl2~2|) . Identical results are obtained with the exact formula 



(26), and are also identical to the results of [11|. Second line corresponds to the calculations 
using harmonic oscillator wave functions with oscillator length a = [(j + I) 1 / 2 /?]" 1 , and 
third line shows the calculations with modified harmonic oscillator wave functions as given 
in Ref.|il. 
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Table 3: Franck- Condon overlap integrals /12 between states in different Morse potentials. 
Parameters in potentials are j\ — 5, R\ — 2.67A, (3\ = 0.90A 1 and j 2 = 5, R 2 = 3.60A, 
f3 2 = 0.60A" . Successive entries correspond to results of formula (^t]) (which agree exactly 
with the calculations of ref.|llj making the integration of the Morse wave functions), the cal- 
culation using harmonic-oscillator wave functions with oscillator length ao = [(J + ^Y^P]^ 1 , 
and the calculation with modified harmonic oscillator wave functions as given in ref.fffl. 
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Table 4: Franck- Condon (squared) overlap between states involved in the transition — 
X 1 S+ of 7 Li 2 molecule. The parameters of this transition are: ji = 69.494, Pi = 0.326 , 
Ri = 5.876 a . ji = 48.106, /3 2 = 0.459 , R 2 = 5.0535 a . Successive entries correspond 
to: Kuschp2[, Drallosf^], Rivas-SilvapOl, Ley-Koo|14| and our results using formula ( p0|) 
with the scaled integral (0). 
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t:.C/UU (J 


i 079 _ i 


(J. -L -L ( -L 












l.OUl o 


1 9^ _ 1 


9 07^ - 1 


1 ^99 _ i 


7 fi07 - 9 


9 Q7^ - 9 




.39406016 


.00180094 


.12553500 


.20752753 


.15219139 


.07606690 


.02974802 


2 


2.103 - 1 


1.511 - 1 


1.345 - 1 


1.000 - 4 


7.110 - 2 


1.374 - 1 


1.263 - 1 




1 OKI _ 1 


8 AAC\ — 9 


1 110 — 1 

1.11U 1 


2Q0 — 2 


9 S9fi — 1 








1.978 - 1 


8.210 - 2 


1.143 - 1 


2 890 - 2 










2.011 - 1 


1.966 - 1 


1.043 - 1 


7.137 - 3 


1.181 - 1 


1.612 - 1 


1.165 - 1 




901 078^ 


1 Qfifi49' : l1 


1 049708^ 


0D71 Q7/lfi 

.UU 1 XO / ^±U 


1 1 snqfizin 


1 «i OA^A 


1 1 fi^il 784 


3 


6.980 - 2 


2.711 - 1 


6.300 - 3 


1.508 - 1 


5.500 - 2 


2.400 - 3 


6.310 - 2 






1 fiQ9 1 


o.uuu 


1 *}Q1 1 


9 800 








9.820 - 2 


1.688 - 1 


1.000 - 4 


1.394 - 1 










5.656 - 2 


2.865 - 1 


3.732 - 2 


1.564 - 1 


1.533 - 2 


3.566 - 2 


1.220 - 1 




.05655825 


.28654540 


.03732449 


.15638530 


.01533330 


.03566019 


.12201313 


4 


1.560 - 2 


1.657 - 1 


1.834 - 1 


3.030 - 2 


6.610 - 2 


1.082 - 1 


2.040 - 2 




4.130 - 2 


1.516 - 1 


6.220 - 2 


4.230 - 2 


8.000 - 2 








9.241 - 3 


1.370 - 1 


2.568 - 1 


5.520 - 4 


1.241 - 1 


6.346 - 2 


1.335 - 3 




.00924096 


.13696229 


.25680238 


.00055202 


.12409529 


.06346440 


.00133720 


5 


2.500 - 3 


5.600 - 2 


2.243 - 1 


6.280 - 2 


1.028 - 1 


3.900 - 3 


8.950 - 2 




8.503 - 4 


3.121 - 2 


2.068 - 1 


1.765 - 1 


3.284 - 2 


6.442 - 2 


9.497 - 2 




.00085027 


.03120572 


.20683663 


.17650291 


.03283566 


.06440957 


.09479757 



Table 5: Franck-Condon (squared) overlap between states involved in the transition B l Tl+ — 
X l Y^ of 7 Li 2 molecule. The parameters of this transition are: ji = 62.05, /?i = 0.354 , 



Ri 

to: 



5.550 a . fa = 48.104, f3 2 = 0.459 , R 2 = 5.0535 a . Successive entries correspond 
3ssel [23] , Drallos[^, PalmapQ|, Ley-Koo[^4]1 and our results using formula (^0|) with 
the scaled integral (^2|). 








1 


2 


3 





.4960 


.3370 


.1220 


.0280 




.4990 


.3200 


.1260 


.0399 




.4940 


.3350 


.1260 


.0349 




.4970 


.3236 


.1279 


.0389 




.5047 


.3295 


.1222 


.0339 




.50469840 


.32944876 


.12216458 


.03391230 


1 


.3470 


.0450 


.2940 


.2150 




.3710 


.0456 


.2400 


.1950 




.3600 


.0426 


.2640 


.2070 




.3694 


.0452 


.2444 


.2006 




.3551 


.0476 


.2638 


.2053 




.35506925 


.04762420 


.26376978 


.20533736 


2 


.1220 


.2940 


.0120 


.1510 




.1110 


.3500 


.0122 


.0989 




.1190 


.3180 


.0137 


.1200 




.1140 


.3423 


.0112 


.1036 




.1149 


.3201 


.0127 


.1157 




.11488206 


.32010981 


.01270236 


.11565949 


3 


.0280 


.2150 


.1510 


.0887 




.0174 


.2230 


.2100 


.0905 




.0234 


.2220 


.1740 


.0925 




.0183 


.2250 


.2031 


.0842 




.0223 


.2215 


.1728 


.0943 




.02224589 


.22154115 


.17282552 


.09426904 



Table 6: Franck- Condon (squared) overlap between states involved in the transition A 2 U 2 — 
X 2 T, + of CN molecule. The parameters of this transition are: j\ = 60.4084, (3\ = 1.2167 , 
Ri = 2.3308 a . j 2 = 73.56 2, /3 2 = 1.22 7 1 , R 2 = 2.2152 a . Successive entries correspond 
to: NichollsQ, McCallum@, Waldenstrom^], Rivas0, Ley-Koo0 and our results 
using formula (|20D with the scaled integral (p2|). 








1 


2 


3 


4 


5 


6 





3.40 - 1 


4.06 - 1 


2.00 - 1 


5.00 - 2 


6.00 - 3 


.00000 


.00000 




4.06 - 1 


4.01 - 1 


1.58 - 1 


3.17 - 2 


3.47 - 3 


2.01 -4 


5.72 -6 




3.38 - 1 


4.06 - 1 


1.97 - 1 


5.00 - 2 


7.00 - 3 


.00000 


.0000 




4.10 - 1 


3.98 - 1 


1.62 - 1 


3.42 - 2 


- 


- 


- 




3.99 - 1 


3.96 - 1 


1.64 - 1 


3.61 - 2 


- 


- 


- 




4.40 - 1 


3.91 - 1 


1.40 - 1 


2.59 - 2 


2.67 - 3 


1.54 - 4 


4.70 - 6 




.44008920 


.39144379 


.13979116 


.02585242 


.00266511 


.00015358 


.00000470 


1 


3.23 - 1 


2.00 - 3 


2.12 - 1 


3.01 - 1 


1.34 - 1 


2.70 - 2 


3.00 - 3 




3.27 - 1 


3.71 - 3 


2.85 - 1 


2.77 - 1 


9.18 - 2 


1.41 - 2 


1.07 - 3 




3.24 - 1 


2.00 - 3 


2.12 - 1 


2.98 - 1 


1.31 - 1 


2.70 - 2 


3.00 - 3 




3.31 - 1 


2.90 - 3 


2.74 - 1 


2.76 - 1 


- 


- 


- 




3.34 - 1 


2.80 - 3 


2.65 - 1 


2.77 - 1 


- 


- 


- 




3.28 - 1 


1.39 - 2 


3.11 - 1 


2.58 - 1 


7.71 - 2 


1.11 - 2 


8.23 - 4 




.32752991 


.01385303 


.31130862 


.25822192 


.07714093 


.01109089 


.00082330 


2 


1.90 - 1 


1.03 - 1 


1.13 - 1 


3.90 - 2 


2.73 - 1 


2.10 - 1 


5.90 - 2 




1.64 - 1 


1.59 - 1 


6.59 - 2 


1.05 - 1 


3.06 - 1 


1.63 - 1 


3.41 - 2 




1.90 - 1 


1.03 - 1 


1.13 - 1 


3.90 - 2 


2.73 - 1 


2.10 - 1 


6.10 - 2 




1.66 - 1 


1.59 - 1 


6.88 - 2 


9.56 - 2 


- 


- 


- 




1.70 - 1 


1.64 - 1 


6.60 - 2 


8.98 - 2 


- 


- 


- 




1.50 - 1 


1.92 - 1 


4.26 - 2 


1.42 - 1 


3.02 - 1 


1.41 - 1 


2.75 - 2 




.15030284 


.19173591 


.04256041 


.14151502 


.30231160 


.14140430 


.02748643 


3 


8.80 - 2 


1.77 - 1 


2.00 - 3 


1.61 - 1 


1.00 - 3 


1.84 - 1 


2.59 - 1 




6.67 - 2 


1.93 - 1 


2.25 - 2 


1.50 - 1 


1.11 - 2 


2.59 - 1 


2.26 - 1 




8.80 - 2 


1.78 - 1 


1.00 - 3 


1.62 - 1 


2.00 - 3 


1.80 - 1 


2.60 - 1 




6.69 - 2 


1.97 - 1 


2.21 - 2 


1.52 - 1 


- 


- 


- 




6.63 - 2 


2.07 - 1 


2.50 - 2 


1.42 - 1 










5.56 - 2 


1.99 - 1 


4.62 - 2 


1.32 - 1 


2.93 - 2 


2.76 - 1 


2.03 - 1 




.05560264 


.19857737 


.04615778 


.13209793 


.02933752 


.27596042 


.20331880 


4 


3.60 - 2 


1.45 - 1 


7.40 - 2 


3.10 - 2 


1.13 - 1 


4.60 - 2 


8.30 - 2 




2.44 - 2 


1.29 - 1 


1.22 - 1 


4.67 - 3 


1.53 - 1 


6.94 - 3 


1.76 - 1 




3.60 - 2 


1.45 - 1 


7.70 - 2 


3.20 - 2 


1.14 - 1 


4.80 - 2 


8.30 - 2 




1.84 - 2 


1.17 - 1 


1.49 - 1 


1.60 - 6 


1.60 - 1 


2.78 - 4 


2.07 - 1 




.01835737 


.11726744 


.14897452 


.00000132 


.16030221 


.00027803 


.20656645 



Table 7: Franck- Condon (squared) overlap between states involved in the transition B 3 U g — 
v4 3 £+ of N2 molecule. The parameters of this transition are: j\ = 58.55, j3\ = 1.301 , 



Ri 

to: 



50.994, j3 2 = 1-280 , R 2 = 2.420 a . Successive entries correspond 
Zare [p6fl , Chang |3l| , Waldenstrom p5| , Palma |p9| , Ley-Koo |14| and our 
results using formula (|20|) with the scaled integral fl22|). 



= 2.287 a . j 2 
Jarmain 



3C 



fl\l>2 





1 


2 


3 


4 


5 


6 


5 


1.40- 2 
8.38- 3 
1.40- 2 


8.90 - 2 
6.57 - 2 
8.60 - 2 


1.26- 1 
1.39- 1 

1.27- 1 


9.00- 3 
4.09- 2 
9.00- 3 


8.60- 2 
4.94- 2 
8.80- 2 


4.40- 2 
1.00- 1 
4.30- 2 


1.07- 1 
5.05- 2 
1.04- 1 




5.69- 3 
.00569285 


5.37 - 2 
.05367159 


1.43- 1 
.14280467 


7.00- 2 
.06997892 


2.74- 2 
.02736416 


1.27- 1 
.12678991 


2.84- 2 
.02832830 


6 


5.00 - 3 
2.80 - 3 
5.00 - 3 


4.20 - 2 
2.92 - 2 
4.40 - 2 


1.19 - 1 
9.94 - 2 
1.13 - 1 


6.90 - 2 
1.03 - 1 
6.90 - 2 


4.00 - 3 
2.04 - 3 
5.00 - 3 


1.04 - 1 
9.29 - 2 
1.06 - 1 


3.00 - 3 
4.02 - 2 
3.00 - 3 




1.71 - 3 
.00170792 


2.14 - 2 
.02138991 


8.95 - 2 
.08952954 


1.25 - 1 
.12450143 


1.45 - 2 
.01447182 


7.41 - 2 
.07414661 


6.89 - 2 
.06915838 



Table 7: 



Continued 



